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Abstract 

We compute the linear responses of the Hubble diagram to small scalar perturbations 
in the Robertson- Walker metric and to small peculiar velocities of emitter and receiver. 
We discuss the monotonicity constraint of the Hubble diagram in the light of these re- 
sponses. 
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1 Introduction 



In an expanding Robert son- Walker universe, the kinematics of general relativity implies 
a one-to-one correspondence between the apparent luminosity of a standard candle at 
rest and its red shift. This correspondence, 'the Hubble diagram', is modified by devia- 
tions from maximal symmetry, "anisotropies" , and by 'peculiar' velocities of the candle 
and its observer. We compute these modifications in linear approximation separately 
for anisotropies and peculiar velocities still in a purely kinematical context. Of course, 
anisotropies and peculiar velocities are intimately related, but this relation presupposes a 
gravitational dynamics, like Einstein's equation and the knowledge of the matter content 
of the universe or like inflation. It also presupposes the use of Boltzmann's equation for 
a gas of candles or ((super-) clusters of) galaxies. 



We assume the kinematics of general relativity: 

• The gravitational field is coded in a time-space metric of signature H , we take 

the velocity of light to be one. 

• Massive and massless, pointlike test particles, subject only to gravity, follow timelike 
and lightlike geodesies. 

• Pointlike clocks, e.g. atomic clocks, are necessarily massive. They move on timelike 
curves and indicate proper time r. 

We add the following cosmological hypotheses: 

• We assume that the metric is Robertson- Walker with small scalar perturbations , 



The scale factor a{t) is a strictly positive function of time only, the perturbations 
b and c are arbitrary functions on time-space, both are much smaller than one in 
absolute value. The separation of the perturbations into anisotropies and inhomo- 
geneities makes no sense for closed universes and by abuse we call the perturbations 
collectively anisotropies. We define the function of one variable by s{x) = sinx 
for the sphere, k = 1, where < x < ^ describes the northern hemisphere. We 
put s{x) = X foi' the Euchdean space, = 0, with < x < oo and six) = sinhx 
for the pseudo-sphere, k = —1, with < x < oo. We take the coordinates x^^^f 
dimensionless and call them 'co-moving position', while the scale factor is measured 
in meters. 

• The test particles are (superclusters of) galaxies and photons. The former are at 
rest (t = T, X, 6^, (y9 = constant) plus a small 'peculiar' velocity. Note that in absence 
of anisotropies and peculiar velocities, the proper time is universal for all these 



2 Our hypotheses 



dr'^ = (1 + 2c) dt^ -a^{l + 2b) [dx^ + d9^ + sin^ Odip^]. 
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timelike geodesies and is taken as time coordinate. Under the same conditions, x 
measures the dimensionless, co-moving, geodesic distance of a position from the 
origin at x = 0. 



3 ChristofFels 



We hst the non-vanishing Christoffel symbols in linear approximation in h and c where 
the underlined terms refer to the symmetric case. We denote by dt ordinary and partial 
derivative with respect to t and similarly for the other coordinates. 

tt 

adta + 2adta{b — c) + a'^dtb, 
\adta + 2adta{b — c) + a'^dtb] , 
sin^ 6 \adta + 2adta{b — c) + a'^dfb] , 
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4 Anisotropics in the Hubble diagram 



The Hubble diagram is a two-dimensional parametric plot. The parameter is the time of 
flight of the photon between the emitting galaxy and receiving one, us today. The two 
observables are the apparent luminosity £ and the spectral deformation z. According to 
our model they are functions of the unobserved time of flight, which is therefore treated 
as parameter and eliminated [21 ■ These calculations are feasible to first order in the 
perturbations of the Robertson- Walker metric. 
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4.1 Trajectories of emitter and receiver 



Our first task is to compute how the trajectories of galaxies are perturbed by the anisotropies 
b and c. In the symmetric case and without pecuhar velocities these geodesies are at rest 
with respect to the co-moving coordinates and of course we take the affine parameter p 
to be proper time t, t = p = t, x = Xe, d = 9e, (f = (fe- We denote by an overdot the 
ordinary derivative of the trajectory with respect to its affine parameter. To first order, 
we have 



With the Hubble rate H := dta/a, the deviations rj. satisfy 



Vt + 
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Ve 



+ 2Hri^ + 
+ 2Hr]0 + 
+ 2H7]^ + 
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To first order, the first equation decouples and we get 

dt 

— = 1 - c (r, Xe, Oe, V^e) 

The other three equations produce peculiar velocities. 



dr 



— = — < exp 



2H{t{f))df 



^'d^ c{t) exp 



(2) 



(3) 
(4) 
(5) 
(6) 



(7) 



2H {t{f))df 



df 



and similarly for the perpendicular components. In this section we will ignore peculiar 
velocities, the emitter is held at rest, only its proper time is affected by the perturbations. 
The perturbed proper time of the receiver is given by a similar formula. 



4.2 Trajectories of photons 

We solve (in first order) the geodesic equation of a photon emitted from a galaxy at time te 
and position Xe, Oe^ V^e and received at time to, today, at our position, which of course we 
take in the center of the universe, x = 0- Fortunately, the singularity of the metric tensor, 
equation (^, in the center is only a coordinate singularity. We need the link between the 
time of flight to ~ te and the geodesic distance x- To zeroth order in the anisotropies h 
and c, the trajectory of the photon is given by t = Ce/a, x = — o-e/^^, 9 = ip = with 
Og := a(te). To first order we write 

i=:ae/a + et, x =: -^e/a^ + e^, 9 =■ ^e, <-P =■ e^- (9) 
The geodesic equation becomes: 

it - — Htt - 2a,He^ + {-) dt{h + c) - 2% d^c + 2 C^) H{h-c) = 0,(10) 
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e^ + 2^He^-2^dtb+ ^ d^{b + c) = 0, 



a 



a" 
as 

as 



-2 "-e 



- c) = 0, 



e^-sin-^^s-'^a^(6-c) = 0. 



To first order, the first two equations decouple and we get the solution 

et + aey = — [b - c), 
a 

o 

where b is the function b evaluated along the zeroth order geodesic: 

Up) ■.= b{t{p),X{p),0e,Ve)- 



(11) 



:i2i 



(13) 



(14) 



(15) 



The desired Unk between the time of flight of the photon and its geodesic distance covered 
is given to first order by: 



dt 



a a 



1 

a 



Let us rewrite this equation in terms of the emission time tg, 



l-{b-c){te,X{Q,ee,^e) ° 



a{te) 



o. , dt 



and integrate 



to 



{b-c){t,Xit),ee,^e 

a{t) 



dt. 



(16) 



(17) 



(18) 



To first order and for a fixed direction {9^, ^Pe) we still have a one-to-one correspondence 
between emission time and geodesic distance. This correspondence is of course direction 
dependent. 



4.3 Spectral deformation 

Now we are ready to compute the spectral deformation of the photon emitted at {tg, Xe, 6'e, ^e) 
with period Te measured by the proper time of the emitter Tg and received at (to, 0, •, •). 
Let us denote by Tq the Doppler-shifted period as measured by the proper time of the 
receiver tq. As the period of the photon is infinitesimal with respect to its time of fiight 
we have 

rto /-to+Todt/dro 

Xe= — 7T ^ / ■ (19) 
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Taylor expanding we obtain 



and the spectral deformation, 



Tedt/dre _ To dt/ drp ^^q^ 



z:=^KfT^ = -[l + bo-K]-l, bo:=b{to,0,;-). (21) 

Note that to first order the spectral deformation is independent of the perturbation c. Note 
also that the spectral deformation now depends on the direction via be ■= b(te, Xe, ^e, V^e)- 

4.4 Apparent luminosity 

We suppose known the absolute luminosity L of the standard candle in Joule per second. 
Our hypothesis about photons flying on geodesies implies that the number of photons is 
constant. The energy E of each photon changes as its frequency 1/T. A unit time interval 
Tg during which a certain number of photons are emitted is measured by the proper time 
Tg. The apparent luminosity i is measured in Joule per second and per square meter. Now 
the unit time interval Tq during which we count the received photons is measured by the 
proper time tq. The relation between the unit time intervals is computed by a formula 
similar to equation (fT^: 

S = ^[l + bo-be]. (22) 

We also need the (orthogonal) detector area determined by a given (infinitesimal) solid 
angle dQ in the direction 9e, ^Pe at time to and at (co- moving) geodesic distance Xe- This 
area is measured by the velocity of light times the proper time tq all squared and is given 
by 

dA = als'^{Xe)[l + 2bo]dn. (23) 

Note that to first order in b and c we may still speak about photons propagating in a 
given solid angle. Finally the apparent luminosity is: 

L dn Eo Te L faeV r , , n 

^=^TT^^ = i — rrr^y — ] l-4&o + 26e. 24 
AndAEeTo 47ra2s2(xe) V«o/ 

Note that as for the spectral deformation, the apparent luminosity to first order does not 
depend on the perturbation c but does depend on the direction. 

4.5 Eliminating the time of flight 

Our last task is the elimination of the unobserved parameter, the time of fiight. To this 
end we differentiate the relation between time of flight and spectral deformation, equation 

Z{Q + 1 = — \l + bo- b{te, X(te), Oe, ipe)] , (25) 
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with respect to te- 

dz ao 

ao 



+ 



a(te 



a{te) 1 + 60 - K^e, X{te) , 9e, ipe) 

o o 
-dtb{te, X{te), 9e, V?e) + d^b{te, X{te), 9e, 'fe)/a{te 



From this and equation (fTTj) we get 



dz dig die 
1 



with 



and integrating 



aoH{z) 



l-bo + c{te{z),X{z),9e,iPe) 

-H{z)-\dt - a{z)-^d^)h{Uz), X{z), ^e, Ve) 



X{z) 



di 

ao Jo H{z) 



X{z) =X{z) + S, 



(26) 



(27) 



:28) 



(29) 



Siz) := -60 X{z) + 



c{te{z),X{z),9e,'fe) 



dz 



ao^^(2;) 
^ {dt - a{z)-^d^)b{Uz),X{z),9,,^, 



aoH{z 



dz. 



4.6 Hubble diagram 

Finally the Hubble diagram is to first order in the scalar perturbations b and c: 



L 



Analiz + l)^s^{X{z)) 



2bo-2-{xi^))6{z) 
s 



(30) 



(31) 



Our unit of time is chosen today and here on earth. Therefore we set cq = 0. Likewise our 
unit of length or more precisely the numerical value of the speed of light is chosen here 
and now and we set bo = 0. This shows that the apparently strongest z-dependence of 

o 

the linear correction to the Hubble stemming from the term —boX is a coordinate artifact. 
The remaining terms are weighted averages of c and a derivative of b along the zeroth 
order path of the photon between the standard candle and us today. In the absence of 
particular conspiracies in the perturbations, the ^-dependence of these terms is weak: 



6{z) 



1 

C(te(^mt), X{Zint), 9^, (fe) [Ot - a{Zint) d^)b{te{Zint) , X{Zint) , 6*6, <^e 



aoH{zi^t) 



aoH{z;ntY 



(32) 



for some intermediate value Zi^t € [0, z]. Indeed, a recent fit to the Hubble diagram jH] up 
to z = 1.8 gives H{z) = Ho{z + 

We conclude that the scalar perturbations produce a Hubble diagram which is a band 
in the z£ plane with more or less constant relative vertical width. This relative width is 
of the same order of magnitude as the perturbations. 



5 Peculiar velocities in the Hubble diagram 

In this section our metric is Robertson- Walker without perturbations, b = c = 0. However 
we admit peculiar velocities of emitter and receiver with respect to the co-moving coordi- 
nates X: 0, ip, or, put more physically, with respect to the cosmic microwave background. 
We compute the changes in the Hubble diagram to first order in the peculiar velocities 
divided by the speed of light, which we have set to one. 



5.1 Trajectories 

We take the line of sight in the direction 9 = 7r/2, ip = 0, and decompose the peculiar 
velocities into parallel and perpendicular components with respect to this direction: Vf> 
Ve\\ + Ve±, Vo = Vo\\ + Vo±. Then we get the initial conditions of the emitter at t = 



i=^=V^, ^ = 0, 0=^, (33) 

the initial conditions of the receiver at t = to 

t= — =Jl + v^, x= — ^ = 0, <^= , 34 

dro V ao aoSo 

and the initial conditions of the go-between ai t = te 

dt 1 
t = — =1, x = -- = 0, = 0. (35) 
dp ae 



The connection between geodesic distance covered by the photon and its time of flight is 

dt 

a(t) ■ 



dt 

Xe= (36) 



5.2 Spectral deformation 

To compute the spectral deformation we have a second photon emitted a period later 
with respect to the proper time Te of the emitter. Therefore this photon will be emitted 
at t = te + \/l + v^Te and at position x = Xe + Ve\\Te/ae- It will be received at t = 
^0 + a/1 + VqTq and at position x = "^^on^o/^o- We therefore have 



Xe + = / , -j-r ■ (37) 
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a/I +^^0 + ^0|| 



Taylor expanding as before yields 

T 

To ao ^1 + 1.2 + 

To first order, the spectral deformation, 

ao 



— (1 - (won -^e||)) - 1, 
tie 

depends only on the difference of the parallel components of the peculiar velocities. 



(3^ 



(39) 



5.3 Apparent luminosity 

We suppose that our candle emits its absolute luminosity L isotropicly in its rest frame. 
When moving with velocity Ve — vo, its emission profile with respect to the receiver is 



dA^ 



L 

4:71 



1 




Ve - 


Vo 
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(1- 




- Vo 
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An (1 - (i;o|| -Ve\\)y 



(40) 



where ip is the angle between the line of sight and Ve — vq. Our convention of orientation 
is such that ip = and Vo\\ — Ve\\ := \ve — "Col cos positive when the emitter moves 
towards the receiver in which case the forward emission is enhanced. The deformed 
emission profile, the first of equations ()40|) . is a special relativistic formula and contains a 
difference of velocities at different points. Its first order approximation, the second part 
of (jlUj) . only contains projections of velocities onto a geodesies and makes sense also in 
general relativity. 

Neglecting Lorentz contractions, which are quadratic in velocity, the detector area 
seen by the first photon is dA = aQS^(xe)df^, while the second photon sees 



dA 



2 2 



Xe + 



VpiiTf, 



Vo\\To 
ao 



dQ 



1 - 2(i;o|| - Ve\\) 



ToS'iXe 



dn. (41) 



The term To/{aoSe) is an atomic period divided by the time of flight and can safely be 
dropped. Note also that we do not have to worry about the angle between the detection 
area and the line of sight which for a moving observer optimizing her efficiency deviates 
from 90° by an amount quadratic in her velocity. Therefore to first order the apparent 
luminosity is: 



div dfi ^t; 

dn dA Ee To 
5.4 Hubble diagram 



An al s'^ixe 



^e||)]. 



(42) 



We eliminate the time of flight as in the preceding section and get the Hubble diagram 
with its linear perturbations coming from the peculiar velocities of emitter and receiver: 



i{z) 



L 



Anal{z + iys^{X{z)) 



I + 2aiXiz))ivo\\ - Ve 



(43) 
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The function <j{x) •= 1 ~ X^'{x)/s{x) vanishes identically for flat universes, k = 0. For 
curved universes, k = ±1, (j{x) is small, 

(44) 

We conclude that peculiar velocities do not perturb the Hubble diagram to first order 
if the universe is flat. For curved universes the linear perturbation is small for small 
redshift and grows with z. 

6 Conclusions 

Under the very general kinematical hypotheses outlined at the beginning, the unperturbed 
Hubble diagram has a monotonicity property jl]. Today, this property is respected by 
supernova data. That might change in a foreseeable future and we already look for 
excuses. The first two that come to mind are fluctuations in the absolute luminosity L 
and absorption by dirt along the line of sight. We find it hard to believe that these two 
effects show a z dependence that will mimic a non-monotonicity in the Hubble diagram. 
We rather expect that they will produce a band in the zi plane (single side band for 
dirt). Two other excuses come to mind next, anisotropics and peculiar velocities. After 
the above calculations, we find it hard to believe, that anisotropics and peculiar velocities 
can account for violations of the monotonicity constraint in the Hubble diagram. 
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